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Abstract 



■ Relativistic non-ideal fluid dynamics is formulated in 3+1 space-time dimensions. The equations 
> . 

. governing dissipative relativistic hydrodynamics are given in terms of the time and the 3-space 

o> ■ 

I quantities which correspond to those familiar from non-relativistic physics. Dissipation is accounted 



for by applying the causal theory of relativistic dissipative fluid dynamics. As a special case 
we consider a fluid without viscous/heat couplings in the causal system of transport/relaxation 
equations. For the study of physical systems we consider pure (l+l)-dimensional expansion in 
^ I planar geometry, (l+l)-dimensional spherically symmetric (fireball) expansion, (l+l)-dimensional 



cylindrically symmetric expansion and a (2+l)-dimensional expansion with cylindrical symmetry 



^ , in the transverse plane [fireharrel expansion). The transport/relaxation equations are given in 

\ terms of the spatial components of the dissipative fluxes, since these are not independent. The 

choice for the independent components is analogous to the non-relativistic equations. 

PACS numbers: 05.70.Ln, 24.10.Nz, 25.75.-q, 47.75. +f 
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I. INTRODUCTION 



The study of nonequilibrium properties of matter produced at relativistic high energy 
nuclear coUisons is important to unravel the underlying interactions between the constituents 
of the system. Nonequilibrium effects are of central importance to the space-time description 
of these collisions. 

In real fluids there are internal processes that result in transport of momentum and 
energy from one fluid element to another on a microscopic level. The momentum transport 
mechanisms give rise to internal frictional forces (viscous forces) that enter directly into 
the equations of motion, and that also produce frictional energy dissipation in the flow. 
The energy transport mechanisms lead to energy conduction from one point in the flow to 
another. 

In this work we use equations of fluid flow that explicitly account for the processes 
described above. We adopt a continuum view here and leave the microscopic kinetic-theory 
view to paper II of this work [1] . In paper II we recover essentially the same set of equations, 
but now with a much clearer understanding of the underlying physics. The alternative 
approach also allows us to evaluate explicitly, for a given molecular model, the transport 
coefficients that are introduced on empirical grounds in the macroscopic equations. For a 
review on nonequilibrium models for relativistic heavy ion collisions see [2]. 

To investigate the space-time evolution of high energy relativistic nuclear collisions in a 
more complete 3-1-1 space-time formalism by means of dissipative fluid dynamics requires 
numerical schemes to solve the fluid dynamical equations of motion and the transport (evo- 
lution or relaxation) equations. The existing numerical schemes for solving relativistic ideal 
fluid dynamics can be used provided the relativistic dissipative fluid dynamics equations, in 
their covariant structure, are cast in a more transparent form. 

The dynamics of relativistic non-ideal fluids is of considerable interest in high energy 
heavy ion collisions [3-13]. In particular the recent application of causal relativistic fluid 
dynamics to heavy ion collisions [4] have proven to be excellent in explaining the spectra of 
particles produced at the Relativistic Heavy Ion CoUider (RHIC) [12, 13]. 

The modeling of relativistic matter produced in these collisions is most conveniently 
performed within the 3+1 formalism, where the 4- vectors, tensors, and the equations of 
motion are decomposed with respect to a general space-time coordinate, allowing to express 
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space-time derivatives in a more transparent way as in non-relativistic mechanics. The 3+1 
representation of the equations for relativistic ideal fluids been discussed by many authors 
(see e.g. [14]), in the context of applications to nuclear coUisions. 

Modeling dissipative processes requires non-equilibrium fluid dynamics and irreversible 
thermodynamics. Standard dissipative fluid dynamics derived from standard irreversible 
thermodynamics was flrst extended from non-relativistic to relativistic fluids by Eckart [15] . 
However, like in its non-relativistic counterpart, in the Eckart theory and a variation thereof 
by Landau and Lifshitz [16], dissipative fluctuations may propagate at an inflnite speed 
(causality problem). In addition, small perturbations of equilibria driven by dissipative pro- 
cesses are unstable [17] and finally, no well-posed initial value problem exists [3]. These prob- 
lems in standard dissipative fiuid dynamics originate from the description of non-equilibrium 
states via the local equilibrium states alone, i.e., it is assumed that local thermodynamic 
equilibrium is established on an infinitely short time-scale. This is a consequence of the 
assumption that the entropy 4-current includes only terms linear in dissipative quantities. 

Causal theories of dissipative fiuids based on Grad's 14-moment method [18] were formu- 
lated by Miiller [19] and generalized relativistically by Israel and Stewart [20]. They were 
formulated to remedy some of the problems encountered in standard theories, in particular 
the causality problem. They are based upon extended irreversible thermodynamics, where 
the entropy 4-current vector of standard thermodynamics is extended by including terms 
that are quadratic in the dissipative quantities. Hence the causal theory is also referred to as 
the second-order theory of dissipative fiuid dynamics. In causal theories of dissipative fiuids 
the set of thermodynamic variables is extended to include the dissipative variables. The 
resulting transport equations and hence the equations of motion arc causal and have well- 
posed initial value problem for realistic equations of state, transport coefficients, relaxation 
and coupling coefficients. Besides the Miiller-Israel-Stewart causal theory of relativistic non- 
ideal fiuids another causal relativistic theory for dissipative fiuids exists, namely, a theory 
by Liu et al [21] which is of divergence type, where the dissipative fiuxes are subject to a 
conservation equation. In this work we will follow the former. 

Causal theories of relativistic dissipative fluids provide the easiest extension of standard 
theory towards flnite signal speeds, and are therefore appealing for nuclear collisions mod- 
eling. Extended fluid theories are still rarely applied in the simulations of nuclear collisions, 
which is also due to the lack of an appropriate formulation. However there are developments 
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in this direction [4, 10]. The purpose of this paper is to provide such a formulation for 
dissipative relativistic fluid dynamics. 

This paper provides a complete set of equations for dissipative fluid dynamics in their 3+1 
dimensional representation, using a causal description of thermodynamics. Having at hand 
an appropriate 3+1 space-time formulation for dissipative relativistic fluid theories, we can 
then employ them in modeling relativistic nuclear collisions. In section II the basic elements 
of causal dissipative hydrodynamics are presented in their covariant compact form. In section 
III we give the causal system of transport equations in their covariant compact form. In 
section IV we give the weak coupling limit of the transport equations and the linearized form 
of the complete system of the 14 equations. In section V wc give the 3+1 representation 
of relativistic dissipative fluid dynamics. As special cases for applications to the study of 
nuclear collision dynamics we specify the system of equations given in Section V to the case 
of (i) (l+l)-dimensional expansion in planar geometry, (ii) (l+l)-dimensional cylindrically 
symmetric expansion and (iii) (2+l)-dimensional expansion with cylindrical symmetry in 
the transverse plane (firebarrel expansion) in Section VI. The (l+l)-dimensional spherically 
symmetric expansion is left for Appendix A. 

Throughout this article we adopt the units h = c = = 1. The signature of the 
metric tensor is always taken to be g^'^ — diag(+l, —1, —1, —1). Upper greek indices are 
contravariant, lower greek indices covariant. The greek indices used in 4-vectors go from 
to 3 {t, X, y, z) and the roman indices used in 3-vectors go from 1 to 3 {x, y, z). The scalar 
product of two 4-vectors a'* , is denoted by Qui, = a^b^. The scalar product of two 

3- vectors is denoted by bold type, namely, a, b, a - b. The notations A^"''^^ = (A"^ + A^")/2 
and A'"'^' = (A"^ — A'^"')/2 denote symmetrization and antisymmetrization respectively. The 

4- derivative is denoted by da = d/dx". Dot product between two tensors will be denoted 
A^'^B^^ while for spatial components of the tensor A • B = A^^B^^ and for the double 
scalar product between tensors A : B = Yliik-^^^^^^ ■ '^^^ notation b (8) b = b^V and 
similarly b (g) (B • a) = 6' Yl,k B^^a!^ denotes the multiphcation of two vectors. 

II. NON-IDEAL RELATIVISTIC FLUID DYNAMICS 

Non-ideal fluids, like ideal fluids, are described by conservation laws for the particle 
current vector and the energy-momentum stress tensor. A thermodynamic equilibrium state 
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of a non-ideal fluid dynamics theory can be characterized by the five independent dynamical 
fields which describe a reversible thermodynamic process in an ideal fiuid dynamics, and one 
often chooses two scalar equilibrium state variables and the three independent components of 
the preferred 4-velocity. In high energy heavy-ion collisions we generally take the net baryon 
number density and energy density as state variables. A general non-equilibrium state or 
irreversible thermodynamic process is to be characterized by 14 independent dynamical 
fields. These are, in addition to the ideal fields, the projections of the stress-energy tensor 
parallel and orthogonal to the preferred velocity (dissipative fiuxes), e.g. the trace and 
trace-free part of the spatial stress tensor, and the spatial heat fiux vector. 

We choose the average 4-velocity u'^ such that the particle flux in the associated rest frame 
vanishes. This is the Eckart frame [15] and is the natural frame in systems where there exists 
some conserved net charge (see [20] for the alternative energy frame description) . The state 
of the fluid is assumed close to a thermodynamic equilibrium state, characterized by the 
local thermodynamic equilibrium scalars such as the equilibrium net charge density ^eq, 
equilibrium energy density £eq the local equilibrium velocity u^^ which in the Eckart frame 
can be chosen such that only the pressure p deviates from the local equilibrium pressure Peq 
by the bulk viscous pressure U = p—peq, whereas n = ricq and e = £eq- Other thermodynamic 
quantities such as the temperature, chemical potential and entropy are obtained from the 
equation of state and thermodynamic relations. 

We consider a fluid that consists of a single component. The variables of concern are the 
net charge 4-current N^^, the energy- momentum stress tensor T'^'^ and the entropy flux S'^. 
The divergence of energy-momentum tensor and of the net charge 4-current vanishes locally. 
That is, the energy-momentum and net charge are conserved locally. However in general 
the divergence of entropy 4-current does not vanish. The second law of thermodynamics 
requires that it be a positive and nondecreasing function. In equilibrium the entropy is 
maximum and the divergence of the entropy 4-current vanishes. Thus 



d^N^ = , 



(2.1) 



0, 



(2.2) 



das'" > , 



(2.3) 
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where 

= nu" , (2.4) 
T^"" = {e + p + n)uV -{p + n)g'"' + 2q^V^ + n"" , (2.5) 

-p(aoq"n-aiq,7r^''^ . (2.6) 

In the local rest frame defined by = (1, 0) the quantities appearing in the decomposed 
tensors take their actual meanings: n = u^j^N'^ is the net charge density, e = u^T'^''u,y is 
the energy density, p + 11 = — | A^i, T^^^ is the local isotropic pressure plus bulk pressure, 
g'* = UiyT'^^A'^ is the heat flow, n'^" = T^'*'") is the shear stress tensor, and s = Uij,S'^ is the 
entropy density. The angular bracket notation, representing the symmetrized spatial and 
traceless part of the tensor, is deflned by A^^""^ = [i (A^ + A^ A^) - | Ai"" A^r] A"^ . 
The space-time derivative decomposes into 9^ = u"D + with m^Vm = 0- ^his space- 
time derivative decomposition D = u'^df, is the convective time derivative and v'' = ^^'^d^ 
is the gradient operator. The projection onto the 3-space A^^" = g^^ — u^u" = A"'^ is 
orthogonal to u", that is, A^'^Ui, — and is the hydrodynamical 4-velocity of the net 
charge and is to be normahzed such that u^u^ — 1 and therefore u^di^Un — 0. Here P = 1/T 
is the inverse temperature. The ai{e, n) and Pi{e, n) in (2.6) are the second order coefficients 
which are expressed in terms of thermodynamic integrals and therefore are given by the 
equation of state. These second order coefficients are presented in more details in paper II 
of this work [1]. 

Prom the entropy 4-current (2.6) one sees that the entropy density and flux are respec- 
tively given by 

s = u^S'' = s{e, n) - (^PqU^ - P.q^q" + /^sTTa^tt"^^ , (2.7) 
^ A'^'^,^^ = fjq" - Pia^nqi" - ai7r^"'q^) . (2.8) 

Note that the entropy density is independent of the ai while the entropy flux is independent 
of the Pi. The negative sign of the non-equilibrium contributions reflects the fact that the 
entropy density is maximum in equilibrium. The thermodynamic coefficients Pi{e,n) > in 
(2.7) model deviations of the physical entropy density from s due to scalar/ vector/tensor 
dissipative contributions to 5'''. The ai{e, n) in (2.8) model contributions due to viscous/heat 
coupling, which do not influence the physical entropy density. 
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III. RELAXATION TRANSPORT EQUATIONS 



As we mentioned before we shall use the Miiller-Israel- Stewart second order phenomeno- 
logical theory for dissipative fluids [19, 20]. Essentially the extended irreversible thermo- 
dynamics theory rests on two hypothesis: (1) The dissipative flows, heat flow and viscous 
pressures, are considered as independent variables. Hence the entropy function depends not 
only on the classical variables, net baryon density and energy density, but on these dissi- 
pative flows as well. (2) At equilibrium, the entropy function is maximum. Moreover, its 
flow depends on all dissipative flows and its production rate is semi-positive definite. As a 
consequence the bulk pressure il, the heat flow q'^, and the traceless shear viscous tensor 
t:^^ obey the evolution equations [20] 

where lo^" — A^"' A^^ d[i3Ua] is the vorticity, O — \/,j,u^ is the expansion scalar, a^'' — 
Xj^i^v^^ — \j^i^vy^ — I A^*^ \/\u^ is the shear tensor and = ii'^ is the acceleration 4-vector. 
The local enthalpy density \s w — e -\- p. Overdot denotes A^p — u^d\Aap. The transport 
coefficients n, (, and r] denote the thermal conductivity, and the bulk and shear viscous 
coefficients respectively. The quantities 

rn = CPo, = kTPi, = 2r](32 , (3.4) 
Inq = C«o, Iqn = reTao, IqTr = uTai, l-j^q = 2r]ai , (3.5) 

are the relaxation times for the bulk pressure Tn, the heat flux and the shear tensor r,r ; 
and the relaxation lengths for coupling between heat flux and bulk pressure {Inq, Iqn) and 
between heat flux and shear tensor [IqT^, l^^g). The CKj and Pi are presented in Ref. [1]. 



kT 



1 

2 



-CO - ^-CTno, ( ) + g^ 



T 2 ^ \kT^ 



-Iq. V. - Iqn V''n + Tgcu'^'q, , 
+l.qV^^q''^+2r^n^^'^u;'^^ , 



(3.1) 



(3.2) 



(3.3) 
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IV. THE WEAK COUPLING LIMIT OF CAUSAL TRANSPORT SYSTEM 



The complexity of the full evolution equations (3.1)-(3.3) makes their apphcations 
tractable only if certain simplifications are made. A particularly simple set of evolution 
equations results from the assumption that there is no viscous/heat coupling (specifically 
tto = Oil = ). We will also drop the term with the one-half factor by assuming that the 
thermodynamic gradients are small. Finally wc make a further assumption that there are no 
coupling between acceleration and dissipative fluxes. Such is the case in systems where the 
gradients in thermodynamic quantities are small and in systems where there is no coupling 
of velocity components. In relativistic heavy ion collisions such cases are, for example, the 
scaling solution [22], effective (1+1) dimensional expansions (expansions in planar geometry, 
in cylindrical/spherical geometry with symmetry - pure radial expansions). The evolution 
equations resulting from (3.1)-(3.3) under the above assumptions are 

n = —(nE-n) , (4.i) 

= -f^-TT''^) , (4.3) 

where 

He = -CO , (4.4) 

q'e - - u>^^ , (4.5) 

< = 2r^a^^ , (4.6) 

are the standard Eckart thermodynamic fluxes. Equations (4.1)-(4.3) are of covariant rela- 
tivistic Maxwell-Cattaneo form. In contrast to the algebraic constraint equations (4.4)-(4.6), 
the evolution equations (4.1)-(4.3) are flrst order partial differential equations, which assure 
that in the local rest frame the viscous bulk/shear stresses and the heat flux relax towards 
their standard limits {JJe , Qe ? '^e'} time-scales ta- The relaxation times ta follow in 
principle from kinetic theory [1]. 

The conservation laws, that is, the net charge conservation equation, the equations of 
motion, and the energy equation together with the evolution equations (4.1)-(4.3) constitute 
a complete system of hyperbolic flrst order partial differential equations for the solution 
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vector of 14 (=1+1+3+1+3+5) dynamical variables {n, e, u^, 11, q*^, n^^}. This system 
represents a 14-Fields Theory for relativistic dissipative fluids which is causal and stable for 
appropriate equation of state, initial conditions, transport and relaxation coefficients [3] . 

The equations of non-ideal fluid dynamics and in particular those of relativistic dissipative 
fluid dynamics are in general rather complicated. It is therefore useful to have a simple means 
for judging both the relative importance of various phenomena that occur in a flow, and 
the flow's qualitative nature. This is most easily done in terms of a set of dimensionless 
numbers which provide a convenient characterization of the dominant physical processes in 
the flow. Flows whose physical properties are such that they produce the same values of 
these numbers can be expected to be qualitatively similar even though the value of any one 
quantity -say velocity or characteristic length/time -may be substantially different from one 
flow to another. Thus the assumption made in this section (the weak coupling limit) must 
be thoroughly checked on a case by case for the problem under consideration. This is the 
subject of current investigation and will be presented somewhere. 

The analysis of the instantaneous dynamics of a relativistic fluid can be summarized by 
the identity 

= a'^M^ + + a;"" + A'^^ , (4.7) 

which is a generalization of the Cauchy-Stokes decomposition theorem. It shows that at 
each space-time point a fluid is accelerated along its proper time axis, and experiences 
shear, rotation (vorticity), and expansion along its local space axes. Explicit expressions for 

A^'^, a^, 9, uo^'^ , and a^^ in terms of the ordinary velocity v and lab-frame space and time 
derivatives are given in section V. Note that from -u^A^*^ = we have u^a^^ = w/'o"^,^ = 
and U/.u'"' = u^'u^^ = and from A^"" = A/^^A^ and A'^'^A^^ = 3 we have A^^a^"" = 
^'^^(^ixv — 0. Prom the shear tensor we define a shear scalar a from 

CJ^uCJ^'' = \(T . (4.8) 

V. RELATIVISTIC DISSIPATIVE FLUID DYNAMICS IN 3+1 FORMULATION 
A. 3+1 space— time formulation 

In order to find a numerical solution to the relativistic dissipative fluid dynamics equa- 
tions, we must first cast the equations into suitable form. The equations are usually written 
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in a more compact covariant form. The equations in this form are not well suited for solv- 
ing on a computer. First it is unclear what information to specify on the boundary of the 
4 dimensional space-time in order to get a well posed problem. Second, this form of the 
equations does not allow us to pose the type of physical questions we often wish to ask. 
Typically we want to be able to specify a physical system at some initial time and find 
the development of this system as time evolves. Because of these two issues, we break the 
covariance of the equations. The covariance is broken in two steps: (i) First we split the 
4 dimensional space-time coordinates into 3 dimensional spatial coordinates plus 1 dimen- 
sional time coordinate, (ii) Secondly, we break a 2nd-rank tensor up into 1-1-3-1-6 (00, Oi, ij) 
independent components. A 4-vector is broken into 1+3 (0,i) components while a scalar 
is unchanged. The numerical problem of constructing the 3-1-1 space-time is (i) Determine 
a set of initial conditions that satisfy conservation equations and transport equations, (ii) 
Evolve this set of initial conditions forward in time. 

In this section we will give explicit expressions for the net charge 4-current N'^, energy- 
momentum stress tensor T'^'^, net charge conservation d^N^ — 0, energy- momentum conser- 
vation d^T'^'^ in terms of the ordinary velocity and lab frame space and time derivatives. In 
addition we also give the explicit expressions for the transport/relaxation equations for the 
dissipative quantities, namely the bulk pressure equation, the heat flux equation and the 
shear stress tensor equation in terms of the ordinary velocity and lab frame space and time 
derivatives. 

The 4-velocity, metric tensor and projection tensor are decomposed as 



(7,u), u = 7V, 7 = (l-vv) 



(5.1) 



9' 




(5.2) 



(7V,7U, -A), A = g-u(8)u 



(5.3) 



where 7 = (1 — v^) 
decomposed as 



and V is the three velocity of the fluid. The space-time gradient is 



D 



(79*, D), D = u-a, d = {a,}. 



(5.4) 
(5.5) 



(7V9t + 7u- d,7u9t - V), V = A-a, 



(5.6) 
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In the irreducible decomposition 

9(8)u = 0-+CU + ^-dA + a u , (5.7) 

the kinematic properties of the fluid, i.e, expansion scalar, acceleration, shear and vorticity 
tensors, are given by 

O = da + Q-u, (5.8) 
a = 7(9tu + Du, (5.9) 
(T ^ ^Q^u> -^SL^u> -^OA, (5.10) 

cu = laul +[a® . (5.11) 

These kinematic quantities can be written in terms of time derivatives and spatial ones 

e = dtj + 'd, (5.12) 



a = -fdtU + BLs, (5.13) 
(T = -{7(u ® dtu^ + ^Adtj - CTs} , (5.14) 
cu = 7'u(?t u' + cUs , (5.15) 



where 



(T, = (9 ® a^ (g) - -^Ai? , (5.16) 

cUs = [5 (g) ul +t a^ (g) , (5.17) 

a, = Du , (5.18) 

1? = a - u . (5.19) 

A complete spht involves 3+1 representations for the fields N^^, T'^^ , and 5''*, together 
with the equations governing their evolution. The representation of T^*^ and 5''^ involves the 
3+1 formulation of dissipative quantities, i.e., iJ, 5° and q = {g*}, as well as 7r°°,7r°' and 
It = TT*-' of dissipative fluxes (77 , g'' , vr^*^). Not all of them are independent. Orthogonality 
to the fluid velocity field, q^u^ = tt^^-u^ = 0, together with = 0, yields three and five 
independent components of heat fiux and shear tensor, respectively. Analogy with non- 
relativistic physics suggests to generally ehminate g° in favor of q, as well as n'^^ and tt'^' in 
favor of 7t, though different choices might be more appropriate in particular applications. 
For discussion on different choices of the dissipative fluxes see [10]. Orthogonality also 
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naturally holds for dissipative forces \y'^\nT — and a'^". The 3+1 representation of the 
orthogonality relations for and n^'^ is 

= 7(gO-q.v)^?° = q-v, (5.20) 
= 7(7r°° - tt'^' . v) =^ tt™ - 7r°' • V = (7t • v) • V, (5.21) 
= 7(7r°' - 7t • v) ^ 7r°^ = 7t • V , (5.22) 

where Eq. (5.22) was used in the last equality of Eq. (5.21). Finally, = becomes 

7r°° - tr(7t) = =^ tr(7t) = (tt • v) • v , (5.23) 

and allows one to eliminate one component of 7t. Substituting q'^ , tt™ and 7r°' in this way, 
it remains to provide expressions for one scalar function, the bulk viscous pressure iJ, three 
components of the heat fiux (i.e, q), and five components of the viscous stress tensor (i.e., 
7t), as in non-relativistic physics. 



B. Conservation Laws 

The local net charge density is subject to the 3+1 representation of Eq. (2.1), and 
the 3+1 representation of Eq. (2.2) yields conservation laws both for the local total energy 
density and the local momentum density 

= dtAf+Q- {Afv} = , (5.24) 
S^T'^o = dtE + d-{{E + V)v + 7(q - (q ■ v)v) + 7t ■ v - ((7t ■ v) ■ v)v} = , (5.25) 
df.T^'' = dtM. + a • {M ® V - T^g + 7(v (g) q - (q • v)(v (g) v)) + 7t - (7t • v) ® v} =(3)^6) 

Thus the charge 4-currcnt vector N^^ and the stress-energy tensor T'^^ are represented by 
m = (TV, N) and T^"" = {E, M, P) where 

J\/=N^ = n7 , (5.27) 

N = n7V = jVv , (5.28) 

^ = rpOO ^ ^ p^^2 _ -p ^ 27q • V + (7t • v) • V , (5.29) 

M = {T%=,,y,, = (£ + V)^\ + 7(q + (q • v)v) + 7t • v , (5.30) 

P = {T%j=,,y,, = (£ + Vh\ 0v-Vg + 27(v 0q>+7t, (5.31) 
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with V — {p + n) the effective pressure. The physical meaning of the components of the 
energy-momentum stress tensor is as follows: T^^ is the total energy density of the fluid, T^* 
is the energy flux density in the ith direction of the flow, T''^ is the momentum density in the 
ith direction, T*-' is the rate of transport of the ith component of the momentum per unit 
volume through a unit area oriented perpendicular to the jth coordinate axis. For a fluid 
at rest every element of area experiences only a force normal to the surface of the element 
and this force is independent of the orientation of the element. If the fluid is ideal, then 
it is non-viscous and will not support tangential stress even when the fluid is in motion. 
The stress acting across a surface in an ideal fluid is thus always normal to the surface. 
A non-ideal fluid will not support tangential stress when it is at rest, but can do so when 
it is in motion. The energy-momentum stress tensor has the viscous stress tensor. The 
pressure tensor T^^ gives the normal component of the surface force (in the j direction) 
acting on a surface element that is oriented perpendicular to the ith coordinate axis. The 
spatial diagonal components of T^'^ are the normal stresses while the off-diagonal spatial 
components are the tangential stresses (or shearing stresses). 
Prom Eqs. (5.29)-(5.31) we can write M. and P as 

M ^ {E + V)v + 7(q - (q • v)v) + Tt • v - ((7t • v) • v)v , (5.32) 
P = M (g) V - 7(v (g) q - (q • v) (v (g) v)) 7t - (7t • v) (g) V , (5.33) 

and the local net charge density and energy density are found from 

e = (M + 7-^q) • V, (5.34) 
n = (1- vv)^/W. (5.35) 

In ideal fluids one notices that M and v are parallel and one can find v from the simple 
expressions [14]. In non-ideal fluids M and v are in general not parallel. Only in special 
cases where the fluid velocity components decouple can one still finds M and v parallel 
and one can still solve for velocity [4]. In the case where there is strong coupling between 
the velocity components one has to solve for v iteratively from Eqs. (5.32) or from from 
Eqs. (5.34) and (5.35). 
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C. Maxwell- Cateneo limit of the transport equations in the 3+1 formalism 



Due to orthogonality it suffices to provide transport equations only for spatial components 
of the dissipative fluxes, i.e. {il, q,7t}. For particular set, Eqs. (4.1)-(4.3), this yields 

-idtU + 7v • ail = —{HE-n), (5.36) 
Tn 
1 

'ydtq + 7v ■ 9q = — (Qe - q) , (5.37) 

jdtn + jv ■ 97t = — (tTb - 7t) . (5.38) 

The complete 3+1 representation of the constraints, Eqs. (4.4)-(4.6) can be written down 
with the help of the orthogonality relations Eqs. (5.20)-(5.23) 

He = -C0, (5.39) 
= -XT {^\dt\nT -V\nT + Si) , (5.40) 
He = 2r]a . (5.41) 

Note that the weak coupling limit of the causal transport equations can be further simplified 
for special cases such as the scaling solution case. In that case the last term (acceleration 
coupled to heat or shear flux) on the right hand side of Eqs. (5.37)and (5.38) vanishes. 
This is also the case in simple (1+1) dimensional problems as we will see in the following 
section (s). 



VI. PHYSICAL PROBLEMS 



In this section we consider special cases with emphasis on different geometries (coordinate 
systems). These cases represent some of the various scenarios one can find in high energy 
nuclear collisions. In this regard the longitudinal direction is taken to be the z axis in all 
scenarios. We will consider: a pure (l+l)-dimensional expansion along the longitudinal z 
axis in planar geometry, a (l+l)-dimensional cylindrically symmetric expansion along the 
radial direction (r axis) and a (2+l)-dimensional expansion along the longitudinal direction 
{z axis) as well as along the transverse direction with cylindrical symmetry (r axis). 

To evaluate the conservation laws d,^N'^ — and d^T'^'^ = in a particular coordinate 
system we calculate N^^ and T^^'^ components in that coordinate system and convert these 
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components into physical components using Eqs. (B2), (B3), (B4) and (B5) from Appendix 
B. Then we apply Eqs. (B17) and (BIO) from Appendix B. 

For (l+l)-dimensional expansion along the 2;-axis in planar geometry derivatives with 
respect to {x, y) must be identically zero by symmetry. So we need to calculate only terms 
containing derivatives in {t,z). Similarly for (l+l)-dimensional spherical symmetric flow 
the terms in [d/dO) and (d/dcj)) vanish identically while for (l+l)-dimensional cylindrical 
symmetric flow the terms in {d/dcf)) and (d/dz) also vanish identically. So we need to 
calculate only terms containing (d/dt) and (d/dr). 

A. (l+l)-dimensional expansion in planar geometry 

In high energy nuclear collisions one often considers (l+l)-dimensional expansion of the 
produced matter along the beam or collision axis taken to be the z-axis. In this case there 
is only one non-vanishing spatial component of the 4-velocity. In relativistic ideal fluid 
dynamics it is sufficient to decompose the net charge 4-current and the energy-momentum 
stress tensor using only one 4- vector, namely the 4-velocity in addition to the metric tensor. 
However in relativistic non-ideal fluid dynamics due to additional dissipative fluxes that are 
orthogonal to the 4-velocity, the decomposition of the net charge 4-current and the energy- 
momentum stress tensor requires additional 4- vectors in order for the decomposed quantities 
to be expressed in terms of the 4-vectors, the metric tensor and the local scalar quantities. 
Since in the local rest frame the components of the dissipative fluxes are spatial the additional 
4-vectors are space-like while the 4-velocity is time-like. In (l-l-l)-dimensional expansion we 
need only one more 4- vector in addition to the 4-velocity. Note that the additional 4-vectors 
are not needed if one wants to keep the dissipative fluxes in their lab frame tensorial form. 
Only if we want to write them in their local rest frame scalar form in the energy-momentum 
stress tensor will we need additional 4-vectors. This is possible for simple cases such as 
( 1 -|- 1 ) -dimensional expansions . 

Prom the line element 

ds'^ = dt"^ - dx^ - dy^ - dz^ (6.1) 
the Minkowski metric is g^^ — diag(l, —1, —1, —1). For a local observer comoving with the 
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fluid, the two 4- vectors are 



u^{t,z,x,y) = (1,0,0,0) , 
m^(i,^,x,y) = (0,1,0,0) , 



(6.2) 
(6.3) 



in the local rest frame: one is the 4- velocity which is time- like, U/^u^ — 1, and the other one 
is space-like, m^m'^ — —1, and it is orthogonal to the 4-velocity, m'^Ufj, — 0. Hence the local 
rest frame is defined by iu'^ = 5f , rh'^ = 5^). Note that in this simple (l+l)-dimensional case 
there is only one independent component of the heat flux, taken to be the z component. In 
the local rest frame the heat flux is = (0, Q^, 0, 0) = = qfh'^ which we simply denote 
by g, the heat flow. In the local rest frame the shear tensor takes the form 

/O \ 

r^^ 



W — 











rVV 



(6.4) 



\0 T« ^ 

where r^^^ represent the shear stress tensor components in the local rest frame of the 
fluid. There is also one independent component of the shear tensor in this simple (1+1)- 
dimensional problem, taken to be the zz component if^^ — r^^ = tt which we simply denote 
it by TT. The tracelessness property of the shear stress tensor then imphes that in Eq. (6.4) 



TT 



yy 



ryy 



~. Note also that the bulk viscous pressure enters the non-ideal 



part of the energy-momentum tensor as —UA^" and in the local rest frame it is — iJA 



where A 
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The net charge 4-current in the local rest frame is 

Ni" = (n, 0, 0, 0) . 
The energy-momentum tensor in the local rest frame is given by 



rpfXV 



(6.5) 





Q 





o\ 


q 


Vi. 














v± 
















(6.6) 



where Vz — P + '7^ + n , V± — p — 7r/2 -|- U. That is, the net charge 4-current and the local 
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rest frame energy-momentum tensor can be decomposed as 

N" = nu" , (6.7) 
f^''' = (£ + V±) u^u"" - V±g^'' + {V, - V±) m'^m'' + 2q m^'^r^ , (6.8) 

The dynamics of the system can be studied by applying a Lorentz boost with Vz in the 
z direction. Thus the net charge 4-current and the energy-momentum stress tensor as 
measured by an observer with velocity Vz with respect to the fluid configuration are given 

by 

= nu^ , (6.9) 
T'^'^ = (£ + Vi:)u^'u'' - V^g"" + {Vz - V±)m^m'' + 2gm(^«'^) , (6.10) 

where u'^ and are given by 

u''{t,z,x,y) = (7, 7^., 0,0) , (6.11) 
m^{t,z,x,y) = {jvz,j,0,0) , (6.12) 

and Vz is the fluid 3-velocity in the z direction and 7 = (1 — v^) Note that g'^w^ = and 
Q = ^/~Q'\i- Note also that = q^/ q. Explicit components of the net charge 4-current 
and the energy-momentum tensor read 

= 771 , (6.13) 

N' = N%z , (6.14) 

rpoo ^ ^^2_p^^2g7V , (6.15) 

= Wj^Vz + g7'(l + vl) , (6.16) 

rpzz ^ y^^2^2 ^ ^ 257^^^ , (6.17) 

with W = e + Vz- Using Eqs. (6.13), (6.15) and (6.16) the local velocity, energy density and 
net charge density can be obtained from 

TVO^^n, £ = r°0-(r'^^ + g)T;„ = (T°o + T'^. + g . (6.19) 
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That is, 



TOO + Vz 
T 



0x2 



rpOO _ 



n 



(6.20) 

(6.21) 
(6.22) 



The net charge and the energy-momentum conservation equations can be written as 



> dtN^ + d,{N%,) = . 

> l-T^' + §;{iT'' + r.)v. + q}^0, 



d^N^" = 
di,T^^ = 



(6.23) 
(6.24) 

(6.25) 



dt dz 

The Maxwell-Cataneo transport equations in simple one-dimensional expansion takes the 
form 



where 



Tqq + q = 
TttTT + tt = — 2 • 2r](T , 



K,T ( ^ + a 



/ = 
/' ^ 7 



■ d d 
■ d d 



/, 
/, 



dt dz 



(6.26) 
(6.27) 
(6.28) 

(6.29) 

(6.30) 

(6.31) 
(6.32) 



The calculational frame components of the shear tensor n^^ are related to the local rest 
frame components r^*^ as follows 



TT 



TT 



7 T = 7 TT , 



TT 



yy 



r-yy 



TT 

2 



(6.33) 

(6.34) 



which can be read directly from Eqs (6.8) and (6.10). Thus due to its simplistic structure, we 
numerically evolve the transverse components (tt^^ or tt^^) since they do not change under 
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the boost and thus they do not involve the Lorentz gamma factors. Thus we numerically 
solve Eq. (6.28). Then we can obtain the longitudinal component tt^^ from which the local 
rest frame quantity r^^ follows from Eq. (6.33). The local rest frame quantity r^^ can also 
be obtained from the tracelessness of n'^'^ which implies that 

B. (2+l)-dimensional axisymmetric expansions 

(2+l)-dimensional axisymmetric expansions are (2+l)-dimensional in Cartesian coordi- 
nates or planar geometry (the flow 4-vectors such as the 4-velocity are functions of time 
and two spatial coordinates) but they are (1+1) -dimensional in one dimensional cylindri- 
cally /spherically symmetric expansions. In this section we consider systems that exhibit 
cylindrical symmetry though the results can be generalized to include spherical (fireball) 
symmetry for radial expansions. In cylindrical coordinates all derivatives with respect to 
the axial and azimuthal directions are zero, and the 4-vectors are functions of time and 
radial coordinates only. We consider the case where the azimuthal and axial components of 
the 4- flow vectors are zero and the only spatial component is the radial component. 

Note that since this is effectively a (H-l) -dimensional problem as a result of symmetry 
considerations there is only one independent component of heat flux taken to be the radial 
component which in the local rest frame is g'' = (0, Q^e^, 0) = (0, ge^, 0) and there are two 
independent components of the shear stress tensor taken to be the zz and 00 components in 
cylindrical polar coordinates. Here = (cos (f), sin (6) . Wc will derive the equations for the 
case with cylindrical symmetry, and then restate important results so that they also apply 
to spherical symmetry and Cartesian coordinates. 

Starting from the planar geometry (where the 4-velocity has two non-vanishing spatial 
components, namely the x and y component) the local rest frame shear stress tensor is 
transformed into the cylindrical coordinate form 

\ 

-T'^'^]cos0sin0 
T'f'^] sin^ (f) + T't"t' 



(6.36) 
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TT' 



pi/ 



/o \ 

T^^ T^2/ 

T^'^ 
yo J 



fo 

[r^'' - T"^"^] cos^ + T"^-^ [t' 

[r''^ -T'^'^]cos0sin0 [r'"'^ 

\0 



The local rest frame net charge 4-current and energy-momentum stress tensor takes the form 

(6.37) 

(6.38) 



m = ( 


n, 0,0,0) , 
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(6.39) 



where 7^^ = p + il + r^'', ^ p + H + t<^<^, V^^p + n + r^' and Vij^p + n + nj. The 
pressure tensor is defined such that 7^'^ is the net rate of transport, per unit area of a surface 
oriented perpendicular to the jth coordinate axis, of the ith component of momentum. 

Since in this {t, x, y, z) coordinate system we have two non- vanishing spatial components 
of 4- velocity there are now more vectors required to construct a complete set of tensors with 
respect to which T'^'^ should be decomposed. The set of four- vector fields that we need are 
u^, P and m'* which in the local rest frame takes the form 



u''{t,x,y,z) = (1,0,0,0) , 
P(i,x,y,z) = (0,0,0,1) , 
m^{t,x,y,z) = (0,6^,0) . 



(6.40) 
(6.41) 
(6.42) 



with the properties UfjU^ — 1, 1^1^ — m^m^ — —1 and it^P = Ui^iri^ — l^rh'^ — 0. One vector 
is time like and the other two are space-like. The space-like vectors are all orthogonal to each 
other and to the time-like vector. The two additional space-like 4-vectors are now required 
since we have flows in two spatial directions. Also unlike in the ideal fluid dynamics we now 
have the coupling of the two non-vanishing velocity components as a result of dissipative 
effects. Thus in the local rest frame the net charge 4-current and the energy momentum 
tensor can be decomposed as 

N" = nui" (6.43) 
f = (£ + V^)u''u'' + V^gi''' + {V, - V^yr + {Vr - V^)m''wy + 2Qru^''m''^ .(6.44) 
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Relativistic dynamics of the viscous heat-conducting fluid is accomphshed by a Lorentz 
boost in the radial direction. For an observer moving with velocity Vr in the radial direction 
with respect to the fluid configuration the net charge 4-current and the energy momentum 
as measured by such an observer are 

N^' = nu^ , (6.45) 
T/^- = (£ + P</,)m^«" + V4>g'"' + {V, - V^)l^r + {Vr - V^)m^'m'' + 2Qu^^'my^ .(6.46) 

where 

u'^{t,x,y,z) = (7,7^;^e^) , (6.47) 
nt,x,y,z) = (7^^r-,0,0,7) , (6.48) 
m''{t,x,y,z) = i^Vr.^er) . (6.49) 

The relationship between the local rest frame and calculational frame components of the net 
charge 4-current and the energy-momentum stress tensor can be read off from Eqs. (6.44) 
and (6.46). The components of the energy-momentum tensor in the projected coordinate 
system are 

yoo = -Vr + IQ'-i^Vr , (6.50) 

T^'- = yv-i'vr + Q''7'(l + vl) , (6.51) 

rj^rr ^ yy^2^2 ^^^^ 2Q'^''vr , (6.52) 

T't"!' = , (6.53) 
T'' = V, , (6.54) 

with W = e + Vr and the usual components are obtained in terms of these as 



rpXX 


= T'^'^cos'^ + T'^'^sin'^ 


(6.55) 


j^yy 


= T'^'^sin^^ + T^^cos^^ 


(6.56) 


jixy 


= [r''''-T'^'^]cos0sin0 , 


(6.57) 


rpOx 


= T°"cos(/), = T°"sin0 , 


(6.58) 


rpxz 


= T''^cos0, T^^ = T^^sin0 . 


(6.59) 



The components of 4-vectors and 2nd-rank energy-momentum stress tensor presented here 
are physical components, and are identical to the components measured with respect to 
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an orthonormal tetrad in a curvilinear, in this case cylindrical, coordinate system. The 
contravariant or covariant components can be obtained from the transformation rules given 
in Appendix B starting from cylindrical coordinates. 

The local velocity, energy density and net charge density can be obtained from 

£ = yoo - (T^r + Qr)v^, = (T^o + rr)vr + Q^ = 771 , (6.60) 

that is, 

rpOr Qr 

= ' (6-61) 

^ = - Too + p. ' (6.62) 
n = (l-^;2)V27vO . (6.63) 

For the space-time derivatives we use the Eq. (B17). The net charge conservation d^j^N^ = 
and the energy-momentum conservation d^j^T^^^ — can be written as 

d^^m = ^ dtN"^ + ^dr{r''N^Vr) = , (6.64) 
d^T^ = =^ 9tr°° + -l9,{r"T'-°} = , (6.65) 
d^T^"^ = =^ atr°'' + ^drir'^T'-'] - '^T^^ = , (6.66) 

where a = 1,2 for cylindrical and spherical geometry respectively. These equations can be 
simplified and written in a better way, suitable for numerical purposes, as follows 

I'V' + lliVVj^-^iVV, (6.67) 
9 j,„ + I {(T" + v. + ,} = ~ (T" + v.-%, (6.68) 

+ 1 {(r»' + g) + = (r* + g) - ^ (v, - r,) . (6.69) 

Note that the structure of these equations is similar to pure (l+l)-dimensional expansion 
except for the geometric terms on the right hand side. These terms results from the trans- 
formation using Eqs. (B17) and (BIO). 

For the transport equations recall that there is only one independent component of the 
heat flux and two of the shear stress tensor. For the heat flux this translates into the 
equation for q, the heat flow. For the shear stress tensor we choose the zz and the 00 
components. This is just for convenience since these choices makes the calculations more 
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tractable. The simplicity comes from the property that these transverse components of the 
shear stress tensor do not change under the boost and hence they do not pick up the Lorentz 
gamma factors which would otherwise render the solution more difficult. In spite of the fact 
that not all of the equations are independent we shall present them all, since depending 
on the problem under consideration, it may be more advantageous using one set instead of 
the other. The weak coupling limit of the causal transport equations with the cylindrical 
symmetry takes the form 



TnU + n = 


-CO, 








-kT + a 















(6.70) 
(6.71) 
(6.72) 
(6.73) 
(6.74) 



where 
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= 9 + a 
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dr 
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-7 
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-a h -& , 
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(6.75) 
(6.76) 
(6.77) 

(6.78) 

(6.79) 

(6.80) 
(6.81) 



The relationship between the calculational frame and the local rest frame components of the 
shear tensor can be read off from Eqs. (6.46) and (6.44). In particular 





(6.82) 




(6.83) 


TT — T , 


(6.84) 


^4"t> — ^4>4> 


(6.85) 


23 





The orthogonality and tracelessness of the shear stress tensor imphes the following relations 
respectively 

if' hi + TT^^ + TT^^ = , (6.86) 

Including the bulk equation we only have three transport equations. Again these equations 
have the same structure as in (l+l)-dimensional case, except for geometrical terms. Note 
that due to symmetry in the problem it can be shown that in Eqs. (3.1)- (3. 3) 

T,u;'-'^g, = T,A'-°A'^^9[;3ii„] = r,A^-9[^ii,]g^ = , (6.88) 
2r/7r"(^<) = . (6.89) 

This is also the case with spherically symmetric expansions. One can deduce that in (1+1)- 
dimensional cylindrically symmetric flow it is not possible to choose a scalar viscous pressure 
as is the case in (1+1) planar geometry and (l+l)-dimensional spherically symmetric flows. 
This is because viscous effects originate from a tensor, which is not, in general, isotropic 
even for (l+l)-dimensional flows. 

In going from the Cartesian coordinate system to the projected coordinate system one 
notices that this is similar to having started with the pure {t, r, 0, z) cylindrical coordinate 
system in which the dynamics of the system are obtained by applying a Lorentz boost in 
the radial direction with 

u''{t,r,(l),z) = (7,7^^r,0,0) , (6.90) 
m^{t,r,(l),z) = (7^^r-,7,0,0) , (6.91) 

as done in Appendix A for the spherical symmetry case using the mixed tensor T^. Note that 
since in this case we have only one non-vanishing spatial component of 4- velocity, we now 
need only one additional spatial 4- vector for energy-momentum stress tensor decomposition. 

C. (3+1)— dimensional axisymmetric expansion: Firebarrel expansion 

We now consider the relativistic dynamics of non-ideal fluids in firebarell geometry. (3+1)- 
dimensional axisymmetric flows are (3+1) -dimensional with respect to Cartesian coordinates 
i.e., the velocity components are functions of time and all three spatial coordinates, but they 
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are only (2+l)-dimensional in cylindrical coordinate system. All derivatives with respect to 
azimuthal direction are zero, and all 4- vector components such as the 4- velocity components 
are functions of only the time, axial and radial coordinates, {t,z,r). In case without swirl, 
the azimuthal velocity component is zero everywhere. As it is much easier to work with 
(2+1) independent variables than (3+1), for axisymmetric flows, it makes sense to work in 
a cylindrical coordinate rather than Cartesian one. 

Let us consider a system undergoing expansion along the longitudinal direction taken to 
be the z-axis and at the same time expanding along the transverse (radial) direction taken 
to be the r-axis. The 4-velocity in cylindrical polar coordinates takes the form 

u'^{t, X, y, z) — (7, ^yVr cos 0, sin 0, ^v^) , (6.92) 

with 

7 = {l-v')-'/\ (6.93) 

y^^yl+y2, (6.94) 

The net charge density, energy density and momentum density components read 

= 771 , (6.95) 

= (£ + 7^)7^ -V + 2^(q'v, + q'vr) + n'' + tt"^"^ + n'' , (6.96) 

= (e + V)-f\, + 7g'(l + vl) + q'lVrV, + n'-'v, + tt^'^i;, , (6.97) 

TO'- = (£ + V)^^Vr + 7g'(l + vl) + q'-fv,Vr + {t:'' + ii'^'^)vr + n'^'v, , (6.98) 

while the other components needed in the conservation equations can be written in term of 
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the above equations 







(6.99) 


rpzz 


= (£ + V)^'^vl + V + 2q'jv, + TT^' , 

\ / 1 A/ J. 1 1 


(6.100) 


rprr 


= (£ + V)-i''vl + V + 2q'-fVr + tt'-^ , 


(6.101) 


rprz 




(6.102) 






(6.103) 


rpXX 
J. 


— l-i. -/. UvJo \ ± 






= [T'''^ - T*^] sin V + 


(6.105) 


rpxy 


= [T'-^ -T-^-^Jcos^sin^ , 


(6.106) 


rpOx 


= T°"cos0, = T"' sin0 , 


(6.107) 


rpxz 


= r''"cos0, T^" = r''"sin0 . 


(6.108) 



Note that we have replaced 7r°° and 7r°* by their spatial components using the orthogonality 
and the tracelessness of the shear tensor tt'*^. In particular 

^00 ^ ^xx ^ ^yy ^ ^zz ^ ^rr ^ + ^zz ^ ^g_^Qg) 

^Qz ^ ^zz^^^^zr^^ ^ {6.110) 
^Or ^ ^^rr ^ ^4>4>jy^ + ^rz^^ ^ (6.111) 

where the first equation for n^^ comes from the tracelessness of n^" (tt^ = 0) and the 
others comes from the orthogonality condition {n^^u,, — 0). Note also that in the projected 
coordinate system one could replace 7r°° by the spatial components using the orthogonality 
relations giving 

TT^O = 7r'"'"v,' + 2n''VrV, + TT^'vl , (6.112) 
from which the tracelessness conditions now becomes 

rr , 4><l> zz 

^ 27r^'VrV, + — = . (6.113) 

For simplicity and in the rest of this section we will now consider the dynamics of a 
relativistic viscous non-heat-conducting system in this firebarrel geometry. For the numerical 
purposes we rewrite the momentum components of the energy-momentum stress tensor in 
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terms of the energy component and the spatial components in terms of the momentum 
components 

= + V- K + n^^v, + n'^Vr , (6.114) 

j.Or ^ (yOO ^ p _ ^zz^^^ ^ ^rz^^ ^ ^g_^^5^) 

rprr ^ yOr^^ ^ ^ ^ _ ^4>4>y2 _ ^rz^^^^ ^ (^q -^-^Q^^ 

rpzz ^ rpOz^^ ^ p ^ ^ZZ^^2 _ ^ZV^^^^ ^ {<0 .111) 

T'' = T^'Vr + i^'^'lil + l^^'v.Vr , (6.118) 

^ T^-v, + 7r'y-fl - K + 7r'^'^]vrV, (6.119) 

In differential form, the (2+l)-dimcnsional conservation equations for net charge density, 

energy density and momentum density written in cylindrical coordinate system read 

= dtN^ + d,N' + ^dr{rN''} , (6.120) 

d^T""^ = dtT"^"^ + a,r°" + ^drirT"^'} = , (6.121) 

d^T^"^ = dtT'^ + d,T'' + ^dr{rT'''} = , (6.122) 

d^T^""- = dtT"-^ + d,T'' + ^drirT'""} - ^T^^ = . (6.123) 

For numerical purposes we write the above set of conservation equations in terms of the 
conserved quantities in the calculational frame, namely the net charge density, the energy 
density and the momentum densities 

dtN^ = -dJN^v,) - -dJrN%r} = , (6.124) 

r 

a^T°° = -9,{T°%,J - 9,{T°%3r} - -T'-^ , (6.125) 

r 

dtT'" = -d,{{T'%, + V + tt^Vt^ - n''v,Vr} - 9,{r^%, + n''-f^^ + n''v,Vr} 

--{T'^'Vr + 7r"77r + ^""v,Vr} , (6.126) 
r 

--{{T'%r + V+ K + tt'^I/t' - ^"'VrV,} , (6.127) 
r 

where 

rpOz rpOr 

Vsz = Vsr^^ . (6.128) 



27 



The local velocities can be obtained from simultaneously solving Eqs. (6.114) and (6.115) 
while the local net charge density and energy density are obtained from 

n^N{l-v^f'\ (6.129) 
£ = - T^^v, - T^Vr , (6.130) 

For the transport equations we need the bulk pressure equation and the equations for the 
non- vanishing shear stress tensor components that are needed in the conservation equations 

(6.131) 
(6.132) 
(6.133) 
(6.134) 
(6.135) 



where 
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Note that 



dtll + [v^dz + Vrdr]n = 
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-n) , 




— (277(7'''' 




dtTT'^'' + [Vzdz + Vrdr]TT^^ = 


— (277(7- 


-TT-) 


dtTT'' + [Vzdz + Vrdr]7r''^ = 




-tt"'^) 









O ^ da + dzi^Vz) + dri^Vr) + ^ , (6.136) 

r 

a- = ^^i^vr) + ^(2 + f ^.')0 , (6.137) 

a- = v'(7^'.) + ^(1 + l'vl)0 , (6.138) 

= I [V^'M + V'ilVr)] + l^\vz0 , (6.139) 

a^t' = -^ + -6> , (6.140) 
r 3 



V. = [dz + I'Vzidt + Vzdz + Vrdr)] = - V% (6-141) 

V, = [dr + l^Vr{dt + + Vrdr)] = - V' ■ (6-142) 

(6.143) 

^xx ^ ^rr ^Qg2 ^ ^ ^<^^ ^-^2 ^ ^ (6.144) 

^yy ^ ^rr ^ ^ ^<^<^ ^Qg2 ^ ^ (6.145) 

V'^ii'' + V^-"^ = V^'lVr - — . (6.146) 

r 
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The local rest frame dissipative fluxes are obtained from the above evolved quantities 

t:^" = A^A^tt"^ = A^A^r"^ . (6.147) 

Note that due to symmetry considerations in our problem we only need to solve the bulk 
pressure equation and the three independent components of the shear stress tensor. In the 
non-relativistic limit the equations for mass density, energy and momentum conservation 
are 

dtp + d,{fyv,) + ^dr{rpvr) = , (6.148) 

dtE + d,{{E + r,)v, + TT^X} + \dr{r[{E + rr)vr + 'K^'^V,]] = , (6.149) 

dtM' + d,{M'v, + V,} + \{r{M'vr + tt^^} = , (6.150) 

dtW + d,{M'v, + Ti''} + -dr{r[M'vr + Vr]} - -P<^ = , (6.151) 

r r 

where 

(6.152) 

(6.153) 
(6.154) 
(6.155) 
(6.156) 
(6.157) 
(6.158) 

The weak coupling limit of the causal transport equations in non-relativistic limit takes 
the form 

rnii + il = -C9-V , (6.159) 
rgq + q=-A9r, (6.160) 
T^TC + 7t = -2770- . (6.161) 

where here the upper dot stands for the material or Lagrangian time derivative (i.e., d/dt — 
d/dt+v-Q). The explicit transport equations for the nonvanishing shear tensor components. 



E 




€ 


1 

- r-1^' 








= PVr , 


Vr 


= p + il + 7^"^ 




^ p + n + n^\ 




= p + n + TT^^ . 
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in the weak coupling limit are 



with 







(6.166) 




= drVr - -9 • V , 


(6.167) 




= ^9 v , 


(6.168) 


ZT 

a 




(6.169) 



For numerical purposes we write the conservation equations, Eqs. (6.148)-(6.151) 



dtM +dz{NVz)+dr{NVr) = 



r 

rpOr 



r 



where 



dtT"' + dz{T''vz + tt'-^} + S^ir^^;, + Vr} 



VSZ = ^^2 = , 

rpOr rj-iQr 

Vsr = "fr = ■ 



6.162) 
6.163) 
6.164) 
6.165) 



as 

6.170) 
6.171) 
6.172) 
6.173) 



6.174) 
6.175) 

6.176) 



In formulating the relativistic equations in this section we decided to keep the lab-frame 
quantities for the dissipative fluxes. Sometimes the local rest frame equations are more com- 
phcated than the lab frame equations. The local rest frame in this case of an effective (2+1)- 
dimensional flow contains Lorentz contraction/dilation factors and higher powers of velocity 
which makes the calculations and computations more difficult. The transport/relaxation 
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equations are readily written in their lab frame and if one wants the corresponding local rest 
frame quantities one makes the necessary Lorentz transformation. In the non-relativistic 
limit the equations are cast in laboratory frame quantities. 

VII. NUMERICAL ASPECTS OF RELATIVISTIC NONIDEAL FLUID DYNAM- 
ICS 

The equations of relativistic non-ideal fluid dynamics are truly formidable; indeed they 
have never been solved for nontrivial problems except in simplified cases such as scaling 
solutions [3]. Analytic solutions of relativistic fluid dynamics are rare and in most cases of 
interest we must use numerical methods [4]. 

The complicated structure of relativistic dissipative fluid dynamics equations makes it 
unattractive to use in modeling relativistic nuclear collisions. Unlike in non-relativistic 
dissipative fluid dynamics we have time derivatives appearing in more than one places. 
This makes computation more expensive. Otherwise the system of equations given here are 
readily solved by existing numerical schemes that are used to solve fluid dynamics problems. 
Knowledge of the equation of state, initial conditions and the transport coefficients is needed 
to completely solve non-ideal fluid dynamics. Once these arc obtained consistently one would 
then be in a position to do a thorough study of the effects of dissipation in nuclear collisions. 

In this section, we discuss basic aspects of numerical solution schemes for relativistic 
fluid dynamics. For the sake of simplicity, we consider the case of one conserved charge only. 
With the deflnitions Eqs. (5.27) - (5.31) and the conservation laws Eqs. (5.24) - (5.26) the 
conservation equations can be solved numerically with any scheme that also solves the non- 
relativistic conservation equations. There is, however, one fundamental difference between 
the non-relativistic equations and the relativistic ones. In order to solve the latter for 
jV, E, M., the net charge density, energy density, and momentum density in the calculational 
frame, one has to know the equation of state p{e, n) and v. The equation of state, however, 
depends on n, e, the net charge density and energy density in the local rest frame of the 
fluid. One therefore has to locally transform from the calculational frame to the rest frame 
of the fluid in order to extract n, e, M from A/", E, v. In the non-relativistic limit, there 
is no difference between n and J\f, or e and E and the equation of state can be employed 
directly in the conservation equations. Also, the momentum density of the fluid is related 
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to the fluid velocity by a simple expression. The transformation between local rest frame 
and calculational frame quantities is described explicitly below. 

A. Mathematical Structure of the equations of Non-ideal Fluid Dynamics 

In section V we formulated the equations of non-ideal fluid dynamics for (1+1)- 
dimensional expansion in planar geometry, (1-1-1) -dimensional cylindrical symmetric expan- 
sion ((l-l-l)-dimensional spherical symmetric expansion is discussed in Appendix A) , and 
(2-|-l)-dimensional expansion with cylindrical symmetry. We now ask how to solve them. In 
this connection it is instructive to count the number of the variables to be determined and 
the number of equations available to determine them. 

In non-ideal relativistic fluid dynamics we must flnd six fluid variables: n, e, p and three 
components of the ordinary velocity v; in addition we must now find nine dissipative vari- 
ables: the bulk pressure U, the three components of heat flow q, and the flve independent 
(non-redundant) components of the shear tensor 7t. These fourteen variables are related 
by fourteen partial differential equations governing the flow of relativistic non-ideal fluid: 
the flve conservation equations which are the net charge conservation equation, the en- 
ergy conservation equation and the three components of momentum conservation equation; 
in addition we now have nine transport/relaxation equations for the dissipative variables, 
namely the bulk pressure equation, the three components of the heat flux equation and the 
flve independent components of the shear tensor equation. In addition we have constitutive 
relations: we choose the equation of state to be p = p{n, e). 

For (1-1-1) -dimensional expansions we need to determine only seven variables: n, e, p, v, 
n, q, 71 from a total of six differential equations and an equation of state. The number of 
variables increases with the number of spatial dimensions. The system of fourteen unknowns 
and fourteen equations (the 14-Field Theory) can be solved for the spatial variations of all 
unknowns as a function of time once we are given a set of initial conditions that specify the 
state and motion of the fluid at a particular time, plus a set of boundary conditions where 
constraints are placed on the flow. 

The techniques for solving these nonlinear equations are under investigation and will be 
presented somewhere. Analytical methods can yield solutions for some simplifled problems, 
for example, (l-l-l)-dimensional flow with scaling solution [3]. But in general, this approach is 
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too restrictive, and we should recourse to numerical methods for a more realistic description. 
One of the effective numerical methods of solving equations of fluid dynamics is to replace 
the original differential equations by a set of flnite difference equations that determine the 
physical properties of the fluid on discrete space and time meshes. 

Given suitable initial and boundary conditions we follow the evolution of the fluid by 
solving this discrete algebraic system at successive time steps. We must make sure that the 
finite difference equations are numerically stable and an efficient scheme must be found for 
handling shocks, which can produce discontinuities in the solution or between mesh points. 
We need numerical techniques, which not only can handle the full nonlinear equations, but 
are also versatile and flexible enough to (1) permit a detailed description of the microphysics 
of the gas; (2) allow for structural complexities in the transient regimes; (3) be generalized 
easily to include various processes of dissipation; (4) account for departures from local 
thermodynamic equilibrium. This is a subject of current investigation and will be presented 
somewhere. 

B. Transformation between Calculation Frame and Fluid Rest Frame 

In principle, the transformation is explicitly given by Eqs. (5.27) - (5.31), i.e., by flnding 
the roots of a set of 14 nonlinear equations (the non-linearity enters through the equation of 
state p(£, n)). In numerical applications, however, this transformation has to be done several 
times in each time step and each cell. It is therefore advisable to reduce the complexity of 
the transformation problem. Boosting to a frame moving with velocity — v (the collective 
velocity of the fluid is v) the calculational and local rest frame transformations for the heat 
flux and shear tensor are given by 



TT' 





(7.2) 
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where 



1/ 



7 7V-^ 
7v I + (7 — l)v~'^vv^ 

r = (0, Q^ Q^ q^) , 

/o 






Vo 



(7.3) 
(7.4) 

(7.5) 



while the bulk pressure iJ transforms like the isotropic pressure p. In Eq. (7.3) v is a 3 
column matrix, v"^ is the transpose of v and I is a 3 x 3 identity matrix. 

First note that unhke in relativistic ideal fluid dynamics M. and v are in general no longer 
parallel. In the ideal case M and v are in parallel and 



M • V = Mv ^ E -e 



where M^IM-IjV^IvI. One then obtains the relationship between e and E 



E-Mv , 



while the relationship between n and J\f 



n 



(7.6) 



(7.7) 



(7.8) 



is the same as for non-ideal fluid dynamics because we are using the Eckart definition of the 

4- velocity. With these equations e and n can be expressed in terms of A/", E, M and v. The 

5- dimensional root search is therefore reduced to finding the modulus of v for given J\f, E, 
and M, which is a simple one-dimensional problem. To solve this, one uses the definition of 
M, 

M^{E + p)v . (7.9) 

This equation can be rewritten as a fixed point equation for v for given J\f, E, M: 

M 

(7.10) 



E + p{E - Mv,MVT^) ' 
The fixed point yields the modulus of the fluid velocity, from which one can reconstruct 
V — V M/M, and flnd e and n via Eq. (7.7). The equation of state p{e,n) then yields the 
flnal unknown variable, the pressure p. 
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Alternatively, one can use Eq. (7.10) to rewrite e and n as 



E- 



n = N 



E + p ' 
1_ 



{E+py 



1/2 



e+p 
E + p 



1/2 



from which we can get e and n by simultaneously solving 

fie,n)^{E-e){E + p)-M'^0, 
g{e,n) = {N^ - n^){E + pf - M^M^ = Q , 



(7.11) 
(7.12) 

(7.13) 
(7.14) 



and an equation of state p{e, n) then yields the final unknown variable, the pressure p. 

However, in the non-ideal case dissipative effects brings complications to the mathemati- 
cal structure of the problem. But in some cases like the (l+l)-dimensional expansions where 
the heat flux and the shear tensor can still be replaced by some scalar quantities, as done in 
section VI, one can still generalize the above results. For example in (1+1) radial expansion 

w - or 



e ^ E 



{E + Vr) ' 



n 



E + Vr' 



1/2 



{E + Vr)\ 

from which we can get e and n by simultaneously solving 



/(£, n)^{E- e){E + Vr) - W - Q 
g{e, n) = {J^ - n^){E + Vrf - J^[M' - QTf = 



(7.15) 
(7.16) 
(7.17) 



(7.18) 
(7.19) 



where Vr — p{s, n) + U + t'^^ and equation of state p{e, n) then yields the final unknown 
variable, the pressure p. In addition one needs the transformations for the heat flux and the 
shear stress tensor which in this case are given by 



7Q^ 



„rr „,2^rr 
TT = 7 T , 



(7.20) 
(7.21) 



where the local rest frame heat flow is QJ" and the local rest frame shear tensor is r^''. 

The complications begin to arise when one moves away from the simple cases like the 
one space dimensional expansions. In two and three space dimensional expansions there are 
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now more than one non-vanishing spatial components of 4-velocity flow. The dissipative 
fluxes, in particular the shear tensor, strongly couples the various components of the non- 
vanishing spatial components of the 4-velocity. Thus making the momentum components 
not to be in parallel with the corresponding non-vanishing spatial components of the 4- 
velocity. In finding the solution to this problem numerically the alternative method discussed 
above for finding the local velocity, energy density and the net charge density is no longer 
plausible. One then uses the generalized first method. In this method the local velocities 
are calculated iteratively from Eqs. (5.32) and then the local energy density and net charge 
density from Eqs. (5.34) and (5.35). In addition the transformation of the heat flux and the 
shear tensor should be calculated from Eqs. (7.1) and (7.2). But the latter transformation is 
not necessary if the energy-momentum stress tensor is written in terms of the calculational 
frame dissipative fluxes since their equations are already in this frame. In more complex 
situations this will be preferred to avoid further difficulties. The calculational frame is then 
the natural frame for the dissipative fiuxes. In the nonrelativistic limit the equations should 
reduce to the lab frame equations of non-ideal Newtonian fiuids as shown in this work (see 
Appendix C). 

C. Numerical Schemes 

To handle the coupled system of the partial differential equations presented here we 
replace the differential equations by suitable discrete approximations and solve these nu- 
merically. In general, to model a heavy-ion collision with non-ideal fluid dynamics requires 
to solve the 5 conservation equations and the 9 transport /relaxation equations in three 
space dimensions. Since this is in general a formidable numerical task, one should resorts to 
the so-called operator sphtting method, i.e., the full 3-dimensional solution is constructed 
by solving sequentially three one-dimensional problems. More explicitly, all conservation 
equations are of the type 



U being jV, E, or M and F{U) being N Eq. (5.27), M Eq. (5.32) or P Eq. (5.33). Such 
an equation is numerically solved on a space-time grid, and time and space derivatives are 
replaced by flnite differences. The solution to the partial differential equation (7.22) in three 
space dimensions is obtained by solving a sequence of one space dimension partial differential 




(7.22) 
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equations 

dtU + diFi{U)^0 , (7.23) 

where i = x. The 3-divergence operator in Eq. (7.22) is spht into a sequence of three partial 
derivative operators. Physically speaking, in a given time step one first propagates the fields 
in X direction, then in y direction, and then in z direction. In actual numerical applications 
it is advisable to permutate the order xyz to minimize systematical errors. The advantage 
of the operator splitting method is that there exists a variety of numerical algorithms which 
solve evolution equations of the type Eq. (7.23) in one space dimension (see, for instance, 
[23] and references therein). 

The transport /relaxation equations on the other hand are of convective type 

dtU+ J2 vAU^FiiU) , (7.24) 

i=x,y,z 

U being il, q, or 7t and F{U) being (il^ - n)/-fTn, (q^ - q)/7'^q or (Tt^ - 7t) /-fr^ in the 
weak coupling limit of the causal transport/relaxation equations. These equations are of 
non-conservative convective type. They should be solved by numerical schemes which solves 
the advective equations of the form 

dtU + v.QU = Sources (7.25) 

rather than the conservative form 

dtU + d{C/v} = Sources (7.26) 

unless the transport equations are recast in the form of Eq. (7.26). The actual numerical 
schemes for future calculations are a subject of current investigation and will be presented 
somewhere. 

VIII. SUMMARY AND DISCUSSION 

We have provided a complete set of equations for rclativistic dissipative hydrodynamics 
in their 3+1 representation. Causality has been accounted for by using the extended causal 
description of thermodynamics, relativistically formulated on a phenomenological level. In 
contrast to the conventionally used compressible Navier-Stokes description of non-ideal hy- 
drodynamics, the equations of extended causal thermodynamics guarantee finite propagation 
speeds of heat and viscous signals and yield stable local thermodynamic equilibria. 
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The extended theories provide a model for evolving dissipative fluids in a causal and 
stable manner. In addition, these theories allow one to describe physical systems which 
cannot be modeled by Navier-Stokes-Fourier, in particular situations where the relaxation 
times are not equal to the microscopic collision time. 

In their simplest form (Maxwell- Cattaneo limit) these transport equations describe the re- 
laxation on finite timescales of dissipative fluxes towards the standard Navier-Stokes-Fourier 
equations, which guarantees finite signal propagation. 

A causality preserving formulation is required whenever the thermodynamic timescale 
becomes comparable to the dynamical timescale and therefore the assumption of local ther- 
modynamic equilibrium is not justified. In many problems of interest the inertia due to the 
dissipative contributions to the energy-momentum stress tensor can be neglected. The 3-1-1 
formulation of the corresponding simplified set of equations is given. 

The five conservation laws for particle number, energy and momentum, and the nine 
evolution equations for the thermodynamic fluxes (i.e., the 14-Fields Theory of Relativis- 
tic Non-ideal Fluid Dynamics) form a hyperbohc system of first order partial differential 
equations tractable by numerical methods. 

Ideally one would like to be able to extract the transport coefficients and associated time 
and length scales from a particular model of an equation of state (interactions among the 
constituents). Then study the effects of dissipative nonequilibrium processes on the space- 
time evolution of the system and on the calculated distributions of particles. And finally 
one would like to compare the predicted distributions with those observed in experiments. 

Finally it is noted that non-ideal fluid dynamics may provide a better way of combin- 
ing hydrodynamic approach and the kinetic (sequential scatering) approach required for a 
realistic description of high enrgy nuclear collisions. 

APPENDIX A: SPHERICALLY SYMMETRIC FLOW: FIREBALL EXPAN- 
SION 

For a one-dimensional, spherical symmetric flow the terms in {d/dd) and [d/dcj)) vanish 
identically. When viewed by an observer moving relative to the new coordinates with proper 
velocity v,- in the radial direction, the physical content of space consists of an anisotropic fluid 
of energy density £, radial pressure Vr, tangential pressure V±i radial heat flux q. Thus when 
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viewed by this comoving observer the covariant net charge 4-current and energy-momentum 
tensor in Minkowski coordinates is 



n, 0,0,0) , 

( e q \ 

q Vr 

V± 

V V± J 



(Al) 



(A2) 



Then a Lorentz transformation leads to 



(A3) 

+ q {m'^u^ + m^u'^) , (A4) 



with = (5^; - M^M^, Peff = \{Vr + 2Px), Vr=P+n + TT,Vj_=p+n- 7r/2, m'^ = 

(7,7^^,0,0) and = (7^^,7,0,0). The heat 4-current and the pressure tensor can be 
written respectively as 



(A5) 
(A6) 



The local rest frame net charge 4-current and the energy-momentum stress tensor are ob- 
tained by boosting back with u'^{t,r,4>,9) = u'^{t,r,4>,9) = (1,0,0,0) and m'^(t,r,(l),9) ~ 
rh'^{t, r, (f), 9) — (0, 1, 0, 0). The nonvanishing components of the net charge and the energy- 
momentum tensor 



' r 



nj, N''" = njVr , 

W7' -Vr + 2q^\r = , 

Wj\ + g7'(l + vj) = = TO'- 



rpr 



_rp4 



W'y^vl + Vr + 2q-f^Vr = T' 



(AT) 
(A8) 
(A9) 
(AlO) 
(All) 
(A12) 



with W = e -\-Vr- The local velocity, energy density and net charge density are obtained 
from 

TVO^^n, £ = r™-(r°'- + g)T;„ T"'' = (T™ + n)^;, + g , (A13) 
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that is 

n = (1 - v2)V27vO . (A16) 

The net charge conservation d^N'^ — and the energy-momentum conservation d^TIf — 
can be written as 

= ^ dtN^ + ^drir^'N^r} = , (A17) 
d^T^ = =^ 9tT™ + -^a^r'^T'"} = , (A18) 
d^T^^ = =^ dtT"^' + ^drir'^T^'} - ^T^ = , (A19) 

where o; = 2 for spherical geometry. For numerical purposes we recast the above equation 
as follows 

a^ATO ^ QiNOyX = -N%r , (A20) 
r 

+ a,{(T°° + Vr)vr + g} = + Vr) - «^ , (A21) 

atT°'- + a,{(TO^ + + Vr] = -a-^{T^' + q) - -{Vr - V±) , (A22) 

r r 

The structure of these equations is similar to the one dimensional case with additional 

geometrical terms. 

For the transport equations recall that there is only one independent component of the 
heat flux and of the shear stress tensor. For the heat flux this translates into the equation for 
g, the heat flow. For the shear stress tensor we choose either the 00 or the 99 components. 
This is just for convenience since these choices makes the calculations more tractable. The 
simplicity comes from the property that these transverse components of the shear stress 
tensor do not change under the boost and hence they do not pick up the velocity or the 
Lorentz gamma factors which would otherwise render the solution more difficult. Thus we 
see that in (1+1) -dimensional spherical symmetric flows all the dissipative fluxes can be 
represented by scalar quantities. In the case in which there is no viscous/heat couplings 
(cti's and Z^'s are zero) the weak coupling limit of the causal transport equations in the 



40 



ce-iii(r„e + cT(^ 



cylindrical and fireball geometry takes the form 
Tnfl + n 



T 



1 



(A23) 
(A24) 
(A25) 



where 



f = 
f ^ 7 



+ -0 
r 3 

d d 
d d 



f , 
/, 



e 



+ a , 

r 



d 



d 



(A26) 

(A27) 

(A28) 
(A29) 
(A30) 



The relationship between the calculational frame and the local rest frame components of the 
shear tensor can be read off from Eqs. (A4) and (A2). In particular 



TT — J T — J TT , 



T 



TT 

2 ' 

TT 

"2 ■ 



(A31) 
(A32) 
(A33) 
(A34) 



Including the bulk equation we only have three transport equations. Again these equations 
have the same structure as in one-dimensional case, except for geometrical terms. 

From the energy-momentum conservation, d^Tj^ = 0, in the calculational frame we have 



u'^d^S +{8 + V,s)0 + d^q'^ = {Vr - V±) 



1 ^ 



a^" + qa^rrii^ , 



(A35) 



+ -A^ 



+qm„a'"' + -Oqm^ = , 



(A36) 
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which are the local energy density conservation equation, Eq. (A35), and the local rest frame 
equations of motion, Eqs. (A36). Contracting the second equation with m^u we get 

m^'d^Vr + {Vr - V±)d^m^' - {e + V±)a^m'' + \eq + vl'd^.q - qm''m''a^, = , (A37) 

where 



a — 



29-2^ 
r 



-(7 + 6> = 3 — . 
2 r 



(A38) 
(A39) 



Explicitly, net charge density, the energy and the radial pressure equations are 
h + nO ^0 , 

e+{e + Ves)0 + -q' = —{Vr - V±)a - 2q- , 
73 7 

p; + 7^ = -(£ + V^^a - {Vr - V±hA - ^ + a] 7g , 

where 



(A40) 

(A41) 

(A42) 



a — 

A = 
/ = 



1 

Vr 



[dt'y + vldr'yvr] , 
7 V r/ 

7 [dt + Vrdr] f , 
7^ [Vrdt + dr]f , 



(A43) 

(A44) 
(A45) 
(A46) 



APPENDIX B: PHYSICAL COMPONENTS AND TRANSFORMATION MA- 
TRICES 

Although the covariant and contravariant components of the vectors (e.g, 4-velocity) 
and tensors (energy- momentum tensor) in curved space coordinates (e.g., cylindrical and 
spherical coordinates) are useful for calculations the normalizations brought in by the metric 
make them inconvenient for physical interpretations. More convenient are components on 
orthonormal tetrads carried by the fluid elements. The hne element is given by 

^ g^^dx'^dx" . (Bl) 
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Since the metric is diagonal we will denote the nonvanishing components by g^^^^^^ and 
respectively g{y){v) for its inverse. For a vector , with the abstract contravariant component 
A^^^ we must associate the physical component 



(B2) 



with no summation implied. For covariant components the physical components are related 
to covariant components by the expression 



(B3) 



with no summation implied. To compute the physical components of a tensor the physical 
components of T^^ in terms of its covariant components are 



Similarly for contravariant components we have 



(B4) 



(B5) 



The parenthetic index has no tensorial significance being merely a label. The physical 
components do not of course transform as tensors, but their transformation law can be 
easily deduced. One notices that the diagonal elements of a mixed second rank tensor are 
the same as the physical components. For mixed tensors the physical components are related 
to the abstract components by 





1/2 

T^M — 
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9(v)(u) 


9{u){u) 







1/2 



or 



T(/x, v) 



9 MM 


1/2 

y -'- av 


9 MM 


g(u)iu) 




9MM 



1/2 



(B6) 



(B7) 



• Spherical symmetry 

The metric and the line element in spherical coordinates read 

g''" = diag(l,-l,-r-^-r-2sin-2 0), = diag(l, -1, -r^ -r^ sin^ 0) , (B8) 

ds^ = dt^ - dr^ - r'^dcf)'' - sin" (B9) 
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and the determinant g — —r^ sin while the transformation matrix for the derivatives 
is 



f d_\ 

dt 

d_ 
dx 



d_ 
dy 

d_ 
\d'z I 



( 











\ 



sin cos Q cos ^ cos Q — sin ^ 



sin sin ^ cos sin ^ cos^ 



y cos — sin j 



d_ 
di 

d_ 
dr 

ld_ 
r dd> 



\ 



d 



(BIO) 



V rsm(f)dd I 

The relations between the abstract and physical components of the 4- velocity in spher- 
ical coordinates read (respectively contravariant and covariant) 



u^^ = Ut, u^^^ =Ur, M*-^-* =u^/r, u^^^ =ug/{rsm(j)) . 

Uq — Ut, Ui—Ur, U2—rU(j,, Uz — T SVO. (f) Uq ■ 



(Bll) 
(B12) 



For a symmetric tensor in spherical coordinates we have, for contravariant components, 
the relations between physical and abstract components 

T^^ = T^<^/r^ T^^ = r^,/(r^sin0), T^' = Tee/{rs\n4>f , (B14) 
with analogous expressions for covariant components. 

Cylindrical symmetry 

The the metric and the line element in cylindrical coordinates read 

g^"" = diag(l, -1, -r-\ -1), g^, = diag(l, -1, -r\ -1) , (B15) 
ds'' = de - dr'' - r''d(t>' - dz' . (B16) 
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and the determinant g — while the transformation of the derivatives is 
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dt 
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dx 
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r d(f) 

d_ 



(BIT) 



The relations between the abstract and physical components of the 4- velocity in cylin- 
drical coordinates are (respectively for contravariant and covariant components) 



— Uf, U^^^ — Ur, 



U 



(2) 



,(3) 



uo = ut, ui = Ur, U2 = ru^, Us = Uz 



(B18) 
(B19) 



For a symmetric tensor in cylindrical coordinates we have (for contravariant compo- 
nents) 



rpQQ rp rpW rp rpVl rp I /TilS rp 

1 — J-tti — -'-rri — -'-rtp/ ' 1 — -'-rz i 



(B20) 
(B21) 



with analogous expressions for covariant components. 



APPENDIX C: NON-IDEAL FLUID DYNAMICS IN THE NON-RELATIVISTIC 
LIMIT 

In Cartesian coordinates {t, x, y, z) the governing equations can be written compactly as 
a single equation as 

dt dx dy dz 
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where 





p 












M^'v, + V, 


u = 


My 


, F = 


Myy^ + TT'^y 




M' 








E 







G = 



pVy 



MyVy + Vy 

M'vy + 

{E + Vy)Vy + 7r^2/^^ + ^zy^^ ^ Qy 



H 



(C2) 



M'^Vz + vr^^ 

Myy^ + TT^y 

M'v, + V, 
_ {E + V,)v, + tt^^t;^ + T^y'Vy + _ 

in which M^ = pvi and Vi = p + 11 + ir'^^Sij. The dissipative fluxes are governed by the 
following relaxational transport equation 



dU dU dU dU 
dt ^ dx ^ dy ^ dz 



where tu — {tu, Tq, t.^} and 



U 



n 








Qy 








^xx 








nyy 








^xy _ ^yx 




T^yz ^ ^zy 




^zx ^ ^xz 





-CO 



(C3) 



dx 

ar 



-2^ - l^) 



-2^ i'-t - -P) 
-V 

-Tj 

'I Idz ^ dxl 



9vy_ . dv^ 
dx dy 

dvz_ , dvy_ 
dy dz 



(C4) 
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in which 



Q^dv^ ^ dvy ^ dv^ 



dx ' dy ' dz ^^^^ 
For one dimensional Cartesian coordinates (t, -z), pure cyUndrical coordinates {t,r,4>,z) 
with — Vz — and pure spherical coordinates {t, r, (p, 6) with V(j, — vq — Q the governing 
equations can be written in one single equation 



-T^ + TTT {r"-^Uv} + r— {r"-^F} + ^ + ^ -H = 

dt r'^-i dr ra~l f)r J Fir r 



dr 



(C6) 



where a = 1, 2, 3 for Cartesian, cylindrical and spherical coordinates respectively and in 
particular for cylindrical coordinates 



(C7) 



The governing equations for the dissipative fluxes can be written in a single compact equation 

1 





p 

















u = 




, F = 




, G = 


V 


, 






E 




VVr + Q' 













dU dU 
dt ^ dr 



[U + F{U)] 



(C8) 



where in cylindrical coordinates 



U 



and in spherical coordinates 



U 



' n ' 










, F = 










' n ' 




Qr 






, F = 











-CO 

.dT 

dr 



-K 



-277 - \0) 
-2^ - \^) 



he 



-co 

or 



(C9) 



-2^ - 1^) 
-2rj - |0) 

-2r, - 10) 



(CIO) 
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while in Cartesian coordinates 



n 




-co 








oz 










(Cll) 






ho 








. ho . 





in which 

= ;:iT|:(-"-V) (C12) 

Note that in the Cartesian coordinates {a = 1) the subscript/superscript r is replaced by 
z and the differentiation with respect to r becomes differentiation with respect to z. Also 
the transverse components of tt'^'^, namely tt^^ and tt^^ do not vanish and they are equal to 
— 7r^^/2. In the spherical {a — 3) case the transverse components are the n'^'^ and the tt^^ 
and each equals In the cylindrical coordinates, a — 2, the transverse components 

of 7^^"' are tt'^'^ and n''. 
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